Velocity Distribution in a Viscous Granular Gas 
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We investigate the velocity relaxation of a viscous one-dimensional granular gas, that is, one in 
which neither energy nor momentum is conserved in a collision. Of interest is the distribution of 
velocities in the gas as it cools, and the time dependence of the relaxation behavior. A Boltzmann 
equation of instantaneous binary collisions leads to a two-peaked distribution with each peak relaxing 
to zero velocity as 1/t while each peak also narrows as l/t. Numerical simulations of grains on a 
line also lead to a double-peaked distribution that narrows as l/t. A Maxwell approximation leads 
to a single-peaked distribution about zero velocity with power-law wings. This distribution narrows 
exponentially. In either case, the relaxing distribution is not of Maxwell-Boltzmann form. 



PACS numbers: 45.70.-n,05.20.Dd,05.70.Ln,83.10.Rs 



I. INTRODUCTION 

Velocity distributions in dilute granular gases are 
generically away from equilibrium because the collision 
processes in such gases are dissipative. Even in the 
dilute limit the velocities of different particles may be 
strongly correlated [|, and therefore the usual descrip- 
tion in terms of single particle distribution functions 
may not be sufficient to determine all the properties of 
the granular gas. Nevertheless, the single particle dis- 
tribution P(r, v, t) contains important information, and 
is particularly interesting because in a granular gas it 
typically deviates from the Maxwell-Boltzmann distribu- 
tion. Quite aside from possible spatial inhomogeneity 
effects such as particle clustering, even the single parti- 
cle velocity distribution P(v, t) in general differs from the 
Maxwell-Boltzmann form. The deviations of these single 
particle distributions from the usual "universal" behavior 
in ordinary gases has been a subject of intense interest 
in recent years [1 These deviations can be 
observed in granular gases that achieve a steady state be- 
cause external forcing balances the dissipative collisions 
among particles @, El > or the y c an be observed in un- 
forced gases as they cool down [H GH [H [H Q [l| . 

Typically, gases equilibrate or achieve a steady state 
via collision processes. Conservation of energy and mo- 
mentum imply Maxwellian velocity statistics in three di- 
mensions |lOl . In lower dimensions the situation is more 
complicated; for example, a gas of hard spheres in one di- 
mension simply retains its initial distribution forever. In 
granular media the situation is in any case more complex 
because energy is certainly not conserved |l7j . Most work 
on relaxation in granular materials, in any dimension, 
has focused on the interesting consequences of energy 
non-conservation. A particularly interesting behavior of 
granular gases induced by energy non-conservation is the 
so-called "inelastic collapse," whereby the energy of the 
gas goes to zero in a finite time ^(| ■ I n the "quasi-elastic 
limit" the inelastic collapse is avoided, and one obtains 
non-trivial asymptotic velocity distributions 0, EI 



with features similar to the ones obtained in this paper. 

In the presence of friction, not only is energy not con- 
served, but momentum is not conserved either, leading 
to interesting new relaxation behavior 0, fTsL figt l20j j . In 
fact, we have recently shown [2(J that in the absence of 
conservation laws random linear mixing can lead to ve- 
locity distributions with algebraic or exponential tails, 
with nontrivial characteristic exponents. In general, con- 
servation laws play a crucial role in the universality of 
the usual velocity distribution properties. 

Our focus here is the effect of viscosity, and conse- 
quently, of momentum non-conservation, on the veloc- 
ity distribution as a one-dimensional dilute granular gas 
cools down. The model consists of N grains on a line (or 
a circle, since we will use periodic boundary conditions). 
The grains move freely except during collisions, governed 
by the Hertz potential, 



V(S k , k+1 ) = %\6\l k+1> 
V(S k , k+ i) =0, 



6 < 0, 



5 > 0. 



(1) 



Here 5k,k+i = Vk+i - yk, Vh is the displacement of gran- 
ule k from its equilibrium position, and a is a prefac- 
tor determined by Young's modulus and Poisson's ratio. 
The exponent n is 5/2 for spheres, it is 2 for cylinders, 
and in general depends on geometry. In this paper we 
only consider cylindrical grains, which leads to consid- 
erable simplification while still capturing important gen- 
eral features of the non-Maxwellian distributions. We 
stress that the one-sided (only repulsion) granular poten- 
tial even with n = 2 is entirely different from a two-sided 
(repulsion and attraction) harmonic potential. 

The displacement of the fe-th granule (k = 1, 2, . . . , N) 
of mass to in the chain from its equilibrium position in a 
frictional medium is governed by the equation of motion 



d 2 y k 

'dr 2 
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-a(Vk - Vk+iY 
-a(y k -i - y k y 



^OiVk - Vk+i) 
^Oivk-i - yk)- 



(2) 
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Here 7 is the friction coefficient and 6{y) is the Hcaviside 
function, 9(y) = f for y > 0, 9(y) = for y < 0, and 
0(0) = 1/2. The Hcaviside function ensures that two 
particles interact only when in contact, that is, only when 
the particles are loaded. Note that for periodic boundary 
conditions yk+N = yk- In terms of the rescaled variables 
and parameters 



yk — C n Xk, 



t — — C n i, 



1 
mvo 



(3) 



where C n = (mwg/a) 1 '", Eq. J2J can be rewritten as 

Xk = -JXk - (x k - x k +i) n ~ 1 0(x k - x k +i) 

+ {x k -i - x k T^ 1 e{xk-i - x k ). (4) 

A dot denotes a derivative with respect to t. The velocity 
Vq is an arbitrary choice in terms of which other velocities 
are now expressed, i.e., it settles the energy scale of the 
system. 

Although the problem might appear relatively simple 
because it is one-dimensional and quasi-linear, the one- 
sidedness of the potential leads to analytic complexities 
even in the dissipationless case 0, 0, [2^J , and even 
greater complexities in the presence of dissipation [Tsl 
Il9| . Our purpose in this paper is to explore analytic 
approximations appropriate for low densities. The low 
density feature of these approximations is implemented 
via the assumption that the collisions are always binary, 
that is, that only two granules at a time are members 
of any collision event, and that at any moment of time 
there is at most one collision. 

We explore two low-density analytic approximations. 
One is based on the Boltzmann equation for binary col- 
lisions, presented in Sec. [H] and elaborated in Sec. 1 1 1 1 1 
and the other is an even simpler model that we call the 
Maxwell model, presented in Sec. IIVI The detailed de- 
scription of a binary collision that must be used as an in- 
put in either of these models is also presented in SeclTTI 
To assess the success of each analytic model, a compari- 
son of the approximate model results with numerical sim- 
ulations of the entire chain is carried out in Sec. [V] In 
that section we also summarize the outcome of our anal- 
ysis. 



Since the problem is one-dimensional, one can keep track 
of the precise conditions under which a collision between 
two particles of given velocities will or will not occur, and 
how these events will change the distribution function. 

The binary collision approximation implies that at any 
one time there is at most one pair of loaded grains, so 
we need to consider in detail the fate of a colliding pair 
from the beginning to the end of a collision 0| . At the 
moment of the start of the collision, which we call t = 0, 
the velocities of the grains are u\ and U2- The collision 
ends at the time r = 2tt/ y^8 — 7 2 , when the grains lose 
contact. It is important to notice that for n = 2 the col- 
lision time is independent of the initial condition. This 
is a feature that makes the cylindrical granule geometry 
much simpler than other shapes. Without loss of general- 
ity we take u\ > 112, since we can always relabel particles 
to make it so. The velocities at the moment of separation 
t are found to be 01 



u i = "i — ^ — U2 — 2 — ' 

/ A* 2 + M , M 2 - A* 
"2 = u i — n — + W2 — o — , 



(0) 



where fi = e _7T / 2 . 

In our further analysis we think of the distance trav- 
eled by the granules during a collision as negligible. For 
small damping 7 this distance is (ui + U2)ir/V&, which 
is to be compared to the mean distance between parti- 
cles. The latter can be made arbitrarily large by lowering 
the density. We also take the collision time between two 
granules as instantaneous. This collision time for small 
damping is r w to be compared with the typical 

mean free time of flight of a particle between collisions. 
With these approximations, the only role played by the 
viscosity is dictated by the collision rule ©■ These as- 
sumptions might conceivably be problematic for the most 
energetic particles that may travel a relatively long dis- 
tance during a collision and a relatively short distance 
between collisions, but explicit analysis of these extreme 
events is difficult and probably not important at suffi- 
ciently low densities. 



II. BOLTZMANN EQUATION AND 
COLLISIONS 

Our analytic starting point in the following sections 
is the Boltzmann equation for binary collisions in a spa- 
tially uniform gas, which describes the rate of change of 
the probability distribution of velocities, P(v,t). If u\ 
and U2 are the initial velocities of a pair of particles just 
before a collision and and u[ and u' 2 their velocities just 
after, then the Boltzmann equation is 

^P(v,t) = Jj dm du 2 \u 1 - u 2 \P(ui,t)P(u 2 ,t) 
x [8(v - u[) + S(v — u' 2 ) — S(v — ui) — 5(v — u 2 )] ■ (5) 



III. BOLTZMANN PROBLEM 

We now start with the Boltzmann equation JSJ in 
which the collision rate depends on the difference in the 
velocities of the colliding pair. To write the explicit 
Boltzmann equation for our system, we must collect all 
the instances that lead to a collision, that is, all the in- 
stances for which one or another of the (5-functions in 
Eq. ||HJ can be satisfied. There arc six possible events 
that contribute to a change in P(v,t). We list the con- 
tribution to dP/dt from each event: 

(a) (v, u) > (u[,u' 2 ); v > u 

Implementing the third (5-function in Eq. JSJ) we 



3 



find 



and define 



at 



P(v,t) 



(a) 



-P(v,t) / du(v-u)P(u,t). (7) 



(b) (u,v) — > (u[,u' 2 ); u>v 

Implementation of the fourth <5-function in Eq. © 

leads to 



Ft P ^ 



(6) 



P(v,t) / du(u - v)P(u,t). (8) 



(c) (ui,u 2 ) — > (v,u' 2 ); ui>u 2 

Here we implement the first (5-function in Eq. @ 

to find 



P(v,t) 



xP 



v/fl ~ t 2v 2fi 2 u 

yfl 2 ~fl Ll 2 -H 



(c) 

2v 



(p 2 + fJ-h 



t)P(u,t)du. (9) 



(d) (ui,U 2 ) — > {v,u' 2 ); ui<u 2 

We again implement the first (5-function in Eq. JSJ : 



xP 



2fi 2 u 



2v 



2v (ft 2 + fi)u 



M 2 - V M 2 - M 



t)P{u,t)du. (10) 



(e) (ui,u 2 ) — >(u[ 1 v); m > u 2 

Now we implement the second (5-function in Eq. 101 : 



xP 



2/At 



2r 



/ M = V/" 2 +M M 2 +M 



2r 



fi 2 + fi M + A* 



,t P(«,()i«. (11) 



(f) («i,« 2 ) — ► (■"'n"); w i < w 2 

Again we implement the second (5-function in 
Eq. ©: 



a\ PM 



xP 



! 2 



2v 



2v 



fi 2 + fl /i 2 + /Lt 



2^ 2 u 



^i 2 + /i fl 2 + fi 



(M 2 - M)« 



, i ) (12) 



The full equation for dP(v, t) j dt is the sum of these six 
contributions. While the collisions of types (a) and (b) 
decrease the probability density P(v,t), collisions (c) 
to (f) increase it. 

We assume a scaling solution of the form, 



(13) 



V u 

X = -T7T, V : 



The Boltzmann equation then is 
dP(v,t) j>{t) 



(14) 



(F(x) +xF'{x)) 



(x - y)F(x)F(y)dy 

1 

(y - x)F{x)F{y)dy 



2x 



xF 



_2f?y_ 

\M 2 -M M 2 -M 
2x (fi 2 + 



[i 2 - [i fi 2 - fi 



2.r 



F(y)dy 



2fi 2 y 



/n* VM 2 ~ M /^ 2 - M 
fi A — fi fi z — ji 



2fi 2 y 



2.r 



xF 



xlv? VM 2 + M M + M 2 
2a; (/i 2 — /i)y 



/i 2 + fi fl 2 + fi 



F(y)dy 



2x 



/i 2 + fi fl + [i 2 



^(^-^fA) F[v)dy . (15) 

Defining z = -f 2 — and recognizing that the 

MM MM 

dummy variable 2 can be relabeled as y after the trans- 
formation, one finds that the third term on the right 
hand side of this equation becomes proportional to the 
fourth, and the fifth to the sixth, so that these pairs can 
be combined to yield the simpler expression 

(F(x)+xF'(x)) 



{x - y)F{x)F{y)dy 

(y - x)F(x)F(y)dy 

r x/ " 2 ( 2x 2fi 2 y 



/^ + i7-oo Vm 2 + m m 2 + m 



xF 



2x 



(m 2 - t*)y 



fi 2 + fi fi 2 + /i 



2fi 2 y 



F{y)dy 



2.x 



/' + 1 Jx/n* Vm +m m +m 



xF 



2,c 



(m 2 - v)y 



fi 2 + fi fi 2 + /i 



F(y)dy. (16) 



Since the only explicit time dependence resides in the 
term, this term must be constant. Integrating this 

term leads to cf>(t) = where c is a constant. There- 

fore, asymptotically, 



<j>{t) ~ f 



(17) 



We have not found an analytic solution to Eq. (|16|) . We 
therefore simulate the equation numerically (next subsec- 
tion) and implement a further approximation that leads 
to an analytic solution that we can compare with the 
numerical results (subsequent subsection). 
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A. Simulation of the Boltzmann equation 

We directly simulate the Boltzmann equation using the 
following algorithm: 

1. We start with N grains whose velocities are inde- 
pendently assigned accordingly to an initial distri- 
bution P{v, 0). 

2. We choose one pair of grains with probability pro- 
portional to the modulus of their relative velocity 
and let them collide, using the collision rulc[S] 

3. Time is incremented by twice the inverse of the 
modulus of the pre-collisional relative velocity. The 
factor of 2 accounts for picked pairs that do not 
collide, since our algorithm forces a collision at each 
step. 

4. We iterate steps ^ to [3] many times and for many 
samples. 

In our simulations we took N = 100 and averaged our 
results over 1000 samples. In Fig. ^wc show the result- 
ing velocity distribution for 7 = 0.9 at different times. 
It is clear that the initial symmetric and single peaked 
distribution develops two distinct peaks as it starts to 
collapse to the ultimate equilibrium distribution, a 8- 
function at v = 0. Moreover, in Fig. [21 we show that 
the width of the distribution and the width of each peak 
both decay as 1/t. In the simulation underlying this fig- 
ure the initial distribution was a double peaked Gaussian, 
chosen because it converges quickly. We find the same 
asymptotic behavior for the initially exponential distri- 
bution, namely, that each peak moves inward as 1/t and 
also shrinks as 1/t. The validity of the scaling solution 
Eq. (|13fl is clearly evident in the scaled rendition of the 
velocity distribution shown in Fig. [3] for different times. 



B. Two-particle model 

While the Boltzmann equation (|16[) does not appear 
amenable to analytic solution, we can formulate a sim- 
pler model that may incorporate the main effects of the 



FIG. 1: Velocity distribution obtained using the simulation 
algorithm detailed in Sec. 1111 Al From left to right and top to 
bottom, the panels correspond to time 0, 8000, 16000, 32000, 
64000 and 128000 in the adimensional units used in this paper. 
The initial distribution is a symmetric exponential. Notice the 
change in the scales as time proceeds. 




FIG. 2: Width of the velocity distribution (circles) and the 
width of each peak (stars). The lines represent the best fit of 
the width vs time for the full distribution (exponent —0.998) 
and for each peak (exponent —1.001). 



system. In this model we have an ensemble of rings of 
size I each containing only two particles. In each ring, 
the dynamics is perfectly deterministic: the two balls will 
keep colliding with each other, back and forth, moving 
toward each other with ever decreasing velocities. For 
each one ring, after n collisions, the velocities, which we 
call u n and v n , are obtained by repeated application of 
the collision rule ©, 



H 2n + (-n) n n 2n - (-n) n 

Un = ~ U H V , 



,2n 



(-M)" , t 1 

n "0 



2 it 



(-M) r 



(18) 



-vo, 



5 




c v t 



FIG. 3: Asymptotic behavior of the scaled velocity distribu- 
tion. Different symbols stand for different times (1000 x 2" 
with n = 0, 1, • • • , 13). The constant c in the figure is ar- 
bitrary and was chosen to facilitate comparison with Fig. |U 




FIG. 4: Dashed lines and dotted lines are respectively veloci- 
ties u n and v„ on a ring with two particles, with ito = — «o = I 
and I = f . The solid curves are the envelopes A(t) (upper) 
and B(t) (lower). 



where uq and vq are the initial velocities. Note that each 
of these velocities alternates from positive to negative 
as the particles move in one direction and then another 
in the ring, as illustrated in Fig. 0] The time that has 
elapsed by the n th collision is obtained as follows. Since 
the grains are in a ring, the distance they have to travel 
between two collisions is I. The travel time is l/\u — 
v\, where u and v are their current velocities between 
collisions. Hence, 



| it - v 



ri-l 

Em- 

fc=0 



1 



1 



(19) 



The sign alternation of the velocities in Eq. (|18f) is 
not important in the effort to understand the long time 
behavior of the distribution. We can approximate the 



velocities by envelope functions A(t) (upper curve) and 
B(t) (lower curve), as illustrated in Fig. 0] 



u n = A(t)u + B(t)v , 
v n = B(t)u +A(t)v . 



(20) 



The initial velocities uq and vq are arbitrary, and have 
only been chosen of equal magnitude for purposes of il- 
lustration in the figure. The envelope functions can be 
found by solving Eq. (|19|) for fi n , substituting this into 
Eq. (|18|l ignoring the minus signs in the (— /j,) n factors, 
and setting t n = t. One finds 



Ait) = - 2 



+ 



B(t) = - 

1 

2 



I wo - vo\ ( I 



I 



--It 



|wo - uo| f 1 



I 



I t 



(21) 



I wo - vol ( 1 



I 



--It 



1 -2 



I wo ~Vo\ (I 



--It 



Hence, if initially the particles had a velocity distribution 
P(i>,0), then the velocity distribution as a function of 
time is 



P(v,t) = 

Asymptotically, 
A(t) = -B(t) 
so that 
P(v,t) = 



Siv - A(t)u - B{t)vo) 



xP(m , 0)P(vq, 0)du dv o . 



I 



2 |it - |* 



(22) 



(23) 




sgn(u - w ) 



xP(u , 0)P(vq, 0)du o dv o . 



(24) 



Therefore, P{u, t) consists of two <5-peaks, one at positive 
and one at negative velocities. If the initial distribution 
is symmetric about zero velocity and is properly normal- 
ized, then 



P(v,t) 



1 



2t 




(25) 



The peaks thus move toward the final velocity u = 
as 1/f. This time dependence is in agreement with the 
Boltzmann equation analysis. 



6 



The two 5-pcaks here reflect the fact that the magni- 
tude of the velocity difference between the colliding par- 
ticles eventually becomes independent of the magnitude 
of the initial velocity difference. Clearly, in the Boltz- 
mann equation analysis this is not quite the case and the 
peaks have a finite width as they converge. However, this 
width decreases in time as 1 /t, approaching the behavior 
of the two-particle ring model asymptotically. 



IV. MAXWELL PROBLEM 

In the Boltzmann equation, the rate of collision of a 
pair of particles depends on their relative velocity. The 
"Maxwell problem" further assumes that the pair of col- 
liding grains is chosen randomly, all pairs colliding with 
the same probability at unit rate. The evolution equa- 
tion J5J for the distribution of velocities is then replaced 
by the even simpler form 



x [S(v 



dui du2P(ui,t)P(u2,t) 



6(v 



-S(v - ui) - S(v - u 2 )] ■ 



(26) 



This problem with a completely general linear colli- 
sion rule (of which Eq. © is a special case) was analyzed 
in pnj . In the language used here, if we assume a scal- 
ing form as in Eq. (|13f) . we arrive at an equation par- 
allel to Eq. JSJ) but now with the time-dependent factor 
— 4>(t) j <j)(t) . Again, we argue that this term must be a 
constant, and integration leads to an exponential decay 
of the width of the distribution instead of the power law 
found in the Boltzmann case, 



(27) 



The decay constant depends on the friction coefficient 7. 
In |20( we described in detail how to evaluate a, and in 
Fig. [5] we show the resulting values of a. Furthermore, 
in [2(| we also analyzed the shape of the distribution. We 
found that (regardless of the initial distribution provided 
that the initial average velocity is zero) F(x) has a single 
peak around x = instead of the double-peaked structure 
that we find for the Boltzmann problem and its two-ring 
simplification, and that it has algebraic tails that decay 
as F(x) ~ x~ v ~ x . The exponent v is also a function of 
7, as shown in Fig. [SJ 



V. NUMERICAL SIMULATIONS AND 
CONCLUSIONS 

The two approximations presented in the previous sec- 
tions lead to very different results for the velocity re- 
laxation behavior of a damped one-dimensional granular 



FIG. 5: Exponential decay parameter a (dashed curve), and 
algebraic decay parameter v for the tails of the distribution 
(solid curve), both as a function of viscosity 7. 



system. One (Boltzmann problem) leads to a velocity 
distribution that is double-peaked, with one at positive 
and one at negative velocities, the peaks moving inward 
as 1 jt toward a zero final velocity. The width of the peaks 
also decreases as 1/t. The second model (Maxwell prob- 
lem) leads to a distribution that is single-peaked, with 
a width that decreases exponentially in time and distri- 
bution tails that decay algebraically. Neither of these 
approximations considers the possible effects of the spa- 
tial distribution of granules. 

In order to assess the validity of these approximations, 
and also to get a sense of the possible effects of the spatial 
distribution of granules, we have carried out numerical 
simulations of a full chain of particles. Our collision rules 
are still as indicated in Eq. ©, and, as before, we assume 
collisions to be instantaneous, but now we actually place 
the granules on a line and keep track of their positions 
so that spatial inhomogencitics can occur if the system 
is so inclined. Results for 10000 granules averaged over 
100 simulations are shown in Fig. The particle density 
is 10~ 3 with the particles initially distributed uniformly. 
The initial distribution is taken to be a symmetric expo- 
nential, and 7 = 0.9. As time proceeds, the initial single 
peak splits into two peaks which move inward and be- 
come narrower, as predicted in our Boltzmann analysis. 
Both the inward motion of the peaks and the width of the 
peaks change as 1/t, as in the Boltzmann approximation 
(see Fig. 0. 

We also tested the scaling hypothesis Eq. EH on our 
simulations. In Fig. |H1 we show that, as in the Boltz- 
mann approximation, the scaling works quite well. More- 
over, in this case, the data can be fitted by a double 
r-distribution 



F(x) 



1 



2b 1 T(b 2 ) 



62-1 



(28) 



Note that this distribution has only two (7-dcpcndent) 
free parameters, b\, which sets the scale, and 62, which 



7 



P. 



- A 

K-^T 


A - 
r--, — 


2 -1 


3 1 2 


1 


1 


1 -0.5 


3 0.5 1 


' 1 A 


K 1 1 



-0.5 -0.25 0.25 0.5 -0.3 -0.15 0.15 0.3 
V 

FIG. 6: Velocity distribution obtained from numerical sim- 
ulations on a line with 10000 granules averaged over 100 re- 
alizations. The particle density is 0.001, the initial spatial 
distribution is uniform, and the initial velocity distribution is 
a symmetric exponential. The damping parameter 7 = 0.9. 
From left to right and top to bottom, time is 1000, 8000, 
16000, 32000, 64000 and 128000 in adimensional units. No- 
tice the change in the scales as time proceeds. 
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t 

FIG. 7: Width of the velocity distribution (circles) and the 
width of each peak (stars). The lines represent the best fit of 
the width vs time for the full distribution (exponent —1.009) 
and for each peak (exponent —1.035). The initial velocity 
distribution was a double Gaussian. 



sets the shape. In Fig. [5] we plot this distribution for 
7 = 0.9 (61 = 0.046 and b 2 = 2.637) together with 
the simulation data obtained with a double Gaussian ini- 
tial velocity distribution. We have tested the double re- 
distribution fit for other values of 7, with equal success. 

Comparing the velocity distribution for the Boltzmann 
approximation (Figs. ^ and [3J and the simulations in a 
line (Figs. E| and 0) we note that the Boltzmann model 
converges to the scaling distribution much more rapidly. 
Initially the simulation results exhibit the same behavior 
as the Boltzmann model, but for longer times the bchav- 




-1 -0.5 0.5 1 



cvt 

FIG. 8: Asymptotic behavior of the scaled velocity distribu- 
tion. Different symbols stand for different times (1000 x 2" 
with n = 7, 8, • • • , 13) and the dashed line is the best fit of 
a double gamma distribution. The constant c in the figure is 
arbitrary and was chosen to facilitate comparison with Fig. 01 



ior of the two distributions for small velocities begin to 
differ. In the Boltzmann case, the probability of finding 
a slow granule is much greater than in the simulation. 
We conclude from these results that 

1. The Maxwell problem is an inadequate represen- 
tation of the velocity relaxation properties of the 
granular chain described in Eq. The predicted 
shape of the distribution does not agree with that 
of the physical chain. 

2. The Boltzmann problem in which we have disre- 
garded spatial dependences captures many of the 
essential features of the velocity relaxation, the 
most important and new feature being the ap- 
pearance of two peaks in the velocity distribution. 
The Boltzmann problem and even the simpler two- 
particle simplification of the problem also capture 
the time dependence of convergence of the two 
peaks into a single one at zero velocity. The slow 
(1/i) convergence is due to the fact that the colli- 
sion rate in the Boltzmann picture slows down as 
the gas cools. By way of contrast, the collision rate 
in the Maxwell problem does not, leading to an ex- 
ponential relaxation in time. 

3. The Boltzmann approximation does not correctly 
capture the late time distribution of the slowest 
particles. This is probably due to spatial correla- 
tions that have been ignored in this approximation 
and that arc currently under investigation |23l | . 

One-dimensional momentum-conserving granular 
gases may exhibit "inelastic collapse," whereby the en- 
ergy of the gas goes to zero in a finite time • Whether 
or not this occurs depends on the number of particles N. 
There is a monotonically increasing relation between N 
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and the coefficient of restitution r for the critical value 
of N c (r) above which collapse occurs, with N c — > oo as 
r — > 1. The "quasi-elastic limit" is the limit r — > 1 and 
TV — * oo, but in such a way as to always remain below 
the collapse threshold. In this limit, a double-peaked 
velocity distribution is also observed [Tfl Hll |. In the 
presence of friction there is no inelastic collapse, and 
we always observe a double-peaked distribution. A 
comparison between those results and ours requires an 
understanding of spatial correlations. This analysis will 
be presented elsewhere [23^ . 

Finally, we point out again that even in our most 
complete simulations we have approximated the collision 
events as instantaneous. While we do not believe this 



to be a perceptible source of error, it would be inter- 
esting (but extremely time consuming) to carry our full 
simulations of the model Eq. (@J) with no further approx- 
imations. 
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